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ABSTRACT 

We present mass models for a sample of 30 high-resolution rotation curves of low surface brightness 
(LSB) galaxies. We fit both pseudo-isothermal (core-dominated) and Gold Dark Matter (CDM) (cusp- 
dominated) halos for a wide variety of assumptions about the stellar mass-to-light ratio. We find that the 
pseudo-isothermal model provides superior fits. GDM fits show systematic deviations from the data, and 
often have a small statistical likelihood of being the appropriate model. The distribution of concentration 
parameters is too broad, and has too low a mean, to be explained by AGDM. This failing becomes more 
severe as increasing allowance is made for stellar mass: NFW fits require uncomfortably low mass-to- 
light ratios. In contrast, the maximum disk procedure does often succeed in predicting the inner shape 
of the rotation curves, but requires uncomfortably large stellar mass-to- light ratios. The data do admit 
reasonable stellar population mass-to-light ratios if halos have cores rather than cusps. 

Subject headings: galaxies: kinematics and dynamics — galaxies: fundamental parameters — dark 
matter 



1. INTRODUCTION 

1.1. LSB Galaxies 

Over the last five years the rotation curves of Low Sur- 
face Brightness (LSB) galaxies and the constraints they 
impose on cosmological theories have received much atten- 
tion in the literature. An LSB galaxy is usually defined 
as a disk galaxy with an extrapolated central disk surface 
brightness > 1 mag arcsec"^ fainter than the typical value 
for "normal" High Surface Brightness (HSB) spiral galax- 
ies (Freeman 1970). Colors, metallicities, gas- fractions, 
and extensive population synthesis modelling all support 
the idea that LSB galaxies are unevolved galaxies with 
low (current and past) star formation rates (e.g. van der 
Hulst et al. 1993; McGaugh & Bothun 1994; McGaugh 
1994; McGaugh & de Blok 1997; de Blok, van der Hulst & 
Bothun 1995; van den Hoek et al. 2000; Bell et al. 2000). 
See Bothun, Impey & McGaugh (1997) for a review. 

The observation that LSB and HSB galaxies follow the 
same TuUy-Fisher (TF) relation requires (in the conven- 
tional picture) that LSB galaxies are dominated by dark 
matter (DM) (Zwaan et al. 1995; Sprayberry et al. 1995; 
Verheijen 1997; de Blok & McGaugh 1996). For reason- 



able stellar mass-to-light ratios T^. low surface brightness 
implies low stellar density. Yet, the extended, low sur- 
face density stellar disks cannot be the major contribu- 
tors to the dynamics in LSB galaxies, as no shift in the 
zero-point of the TF relation with surface brightness is 
observed. This contrasts with the dominance of the stel- 
lar population in HSB galaxies of similar luminosity. 

The modest T^, values as implied by the blue colors and 
the (baryonic) TF relation, together with the diffuseness 
of the stellar disks make analyses of the DM distribution in 
LSB galaxies less ambiguous than in HSB galaxies where 
the stellar component can be significant even for fairly low 
T^ values. LSB galaxies are therefore ideal laboratories for 
measuring the distribution of dark matter for comparison 
with predictions of theories of galaxy formation. 

For example, one of the results of numerical Cold Dark 
Matter (CDM) simulations is a so-called "universal halo 
mass density profile" (Navarro, Frenk & White 1996), com- 
monly known as a "NFW profile." NFW (and all GDM) 
mass density profiles are characterized by steep central 
cusps. This is in contrast with the other commonly used 
"classic" pseudo-isothermal sphere halo model which is 
characterized by a constant density core. The parame- 
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ters of the NFW mass density distribution are related to 
the mass of the halo and the density of the universe at the 
time of collapse and are therefore set by the cosmology. 
As these parameters can be determined from observations, 
this opens the possibility of testing the NFW CDM model 
as well its underlying assumptions. 

A first analysis of LSB galaxy Hi rotation curves by de 
Blok & McGaugh (1997) indicated that they did not rise 
as steeply as their HSB counterparts of similar luminosity, 
contrary to CDM predictions. The mass distribution in 
LSB galaxies is more extended and of lower density than 
in HSB galaxies (de Blok & McGaugh 1996). 

Other results also indicate that the steep rotation curves 
implied by CDM are hard to reconcile with the observed 
shallow rotation curves of dwarf galaxies (Moore 1994; Flo- 
res & Primack 1994; Blais-Ouellette, Amram & Carignan 
2001; Cote, Carignan & Freeman 2000). To explain this 
discrepancy the possibility of redistribution of the (cuspy) 
DM due to violent star-formation (thus creating the ob- 
served cores) was sometimes raised, but this has been 
shown to be inconsistent with other observational data 
(Mac Low & Ferrara 1998). 

McGaugh & de Blok (1998) argued that the shapes of ro- 
tation curves of LSB galaxies were inconsistent with those 
predicted by the NFW prescription. This could not be 
explained by the effects of star formation as the larger 
masses of LSB galaxies would require large bursts in or- 
der to redistribute matter on large scales. Their quiescent 
evolutionary history argues strongly against this (van den 
Hoek et al. 2000). 

This comparison with the CDM model is often dismissed 
because of the limited resolution of the observed Hi curves. 
The early Hi LSB rotation curves were obtained using 
the VLA and WSRT radio synthesis telescopes. The rela- 
tively large beams of these instruments resulted in rotation 
curves with only a limited resolution, de Blok & McGaugh 
(1997) did however show that for the best resolved cases 
the effects of beam smearing were not strong enough to 
explain the observed shallow curve as simply the result 
of a steep NFW model curve affected by beam-smearing. 
Similar results were found for more fashionable cosmolo- 
gies, such as ACDM (fi^ ~ 0.3, Oa ^ 0.7), though with 
smaller discrepancies. 

Even so, the theoretical debate now seems to have set- 
tled on halos with cusps even steeper than NFW halos 
(Moore et al. 1999), thus worsening the possible conflict 
between the data and the simulations. From the observa- 
tional point of view the easiest and least ambigous way to 
test the reality of these discrepancies is to measure high- 
resolution rotation curves. 

1.2. Optical rotation curves 

Optical Ha rotation curves of five LSB galaxies from 
the sample of dc Blok, McGaugh & van der Hulst (1996) 
(BMH) were presented in Swaters, Madore & Trewhella 
(2000) (SMT). Though SMT found that for two of the five 
galaxies the inner slopes of the rotation curves were steeper 
than derived from the Hi observations, this difference does 
not affect the BMH conclusion that LSB rotation curves 
have shallower slopes than HSB rotation curves of simi- 
lar amplitude. Because of these steeper slopes, SMT de- 
rive higher maximum-disk T,^ values (in some cases > 10), 



strengthening one of the conclusions from de Blok & Mc- 
Gaugh (1997) that the maximum T^ values in LSB galax- 
ies are too large to be accommodated by reasonable star 
formation histories and Initial Mass Functions. Such high 
values are inconsistent with the existence of a baryonic 
TF-relation (McGaugh et al. 2000). 

A different approach was taken by van den Bosch et al. 
(2000). They attempted to apply a rigorous correction for 
beamsmearing to the BMH Hi data, and thus to derive the 
true "infinite resolution" rotation curve. They conclude 
that the data are not of high enough resolution to accept 
or reject the NFW hypothesis with any significance. How- 
ever, as they use a modified NFW profile with the inner 
slope of the mass-density distribution as an (additional) 
free parameter, it is not clear how significant this conclu- 
sion is. The usual 3-parameter rotation curve fits are al- 
ready under-constrained; adding another parameter does 
not improve the significance of the results. Furthermore, 
in some cases they find such low values for the inner slope 
that their NFW-halos effectively become core-dominated. 
These halos do of course fit the data, but do not occur in 
CDM simulations. 

The general picture as derived from early observations 
of rotation curves of LSB galaxies therefore still holds: 
LSB galaxies are unevolved, low density galaxies, domi- 
nated by DM. Their rotation curves have shallower slopes 
than those of HSB galaxies of similar amplitude, and the 
shape of the best-resolved LSB curves are not necessarily 
consistent with the NFW rotation curve shapes. 

1.3. New data 

In this paper we present an analysis of high-resolution 
high-quality hybrid Ha/Hi rotation curves of a sample of 
30 LSB galaxies. Of this sample 26 curves were taken 
from the large sample of 50 LSB galaxies presented in Mc- 
Gaugh, Rubin & de Blok (2001) (herafter Paper I). In 
that paper an extensive description is given of the data, 
the sample and reduction method. We also refer to that 
paper for a comparison of the new Ha data with the BMH 
Hi curves. We also re-analyse the data for an additional 5 
curves taken from Swaters, Madore & Trewhella (2000). In 
this paper we derive mass models under various assump- 
tions for Ti, and fit these models both with NFW halos 
and pseudo-isothermal halos. A similar analysis for a dif- 
ferent set of rotation curves of dwarf and LSB galaxies is 
given in de Blok & Bosma (2001). 

In Section 2 we discuss the sample, and discuss the 
derivation of the rotation curves. We also show internal 
and external comparisons of the data and discuss possi- 
ble systematics. In Section 3 we discuss the various mass 
models. Section 4 contains the results of the model fit- 
ting. Section 5 discusses the implications for the various 
halo models. In Section 6 we turn our attention to the 
maximum disk and a summary is given in Section 7. 

When using absolute distances we have used a Hubble 
constant Hq — 75 km s^^ Mpc^^. 

2. THE DATA 

2.1. Sample and raw data 

The data and reduction methods have been extensively 
described in Paper I. In summary, we use long-slit major 
axis spectra taken with the 4-m telescope at Kitt Peak in 
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June 1999 and February 2000 and the 100" telescope at 
Las Campanas in November 1998. Velocities were derived 
from the intensity weighted centroid of the Ha and [Nil] 
lines. 

As the aim of this excercise is to derive mass models 
which can yield significant constraints on the distribution 
of DM we select only the 26 high quality galaxies from Pa- 
per I. We split this high-quality sample in two sub-samples. 
Sample I contains LSB galaxies from BMH and van der 
Hulst et al. (1993) for which a full set of photometry and 
Hi data is available. In Sample I we also include the five 
galaxies presented in SMT. For these galaxies Hi and op- 
tical photometry are taken from BMH. Tables 1 and 2 
contain a full list of the galaxies analysed here, along with 
some of their global parameters. 

Sample H consists of ESO-LV and UGC LSB galaxies, 
for which an optical rotation curve is available but no op- 
tical or Hi photometry. 

For the galaxies in Sample I, Hi observations are avail- 
able which often extend to larger radii than the Ha data. 
To make the best use of both types of data, we have con- 
structed hybrid rotation curves. These consist of the Ha 
data over the range of radii where available, and 21-cm 
data to define the outermost points. No attempt has been 
made to "average" the different types of data: Ha is given 
precedence over the range of radii where it is available. 

2.2. Derivation of the smooth curves 

One of the main assumptions made when deriving mass 
models from rotation curves is that the gas and stars trace 
circular orbits in an axisymmetric potential. Though the 
shape of the optical rotation curves in Paper I is well- 
defined, the scatter between individual datapoints means 
we cannot simply use the raw rotation curves to estimate 
the radial run of the gravitational potential. For this one 
needs a smooth curve which retains real small-scale details, 
but without the observational scatter. 

The method most often used to produce these smooth 
curves is to fit splines to the data. Here we use have used 
a robust version of this procedure (local regression, see 
Loader 1999). The smooth curves were rebinned to a bin 
width of 2". The error bars in the rebinned data points 
consist of two components: one due to observational errors 
caused by the measurement uncertainties in the individual 
raw datapoints (for this we use the average weighted mea- 
surement error in each bin), and an additional component 
caused by differences between approaching and receding 
sides and non-circular motions (which we define as the 
difference between the weighted mean raw velocity and 
the velocity implied by the spline fit at that radius). For 
the final error estimate these two uncertainties were added 
quadratically 

For some high S/N data points the error bars become 
unrealistically small (sometimes less than 1 km s^^). This 
has no physical significance and simply tells us that the 
profile centroids were well-determined. These small error 
bars can however easily dominate any model fit and can 
severely bias x^ values or goodness-of-fit parameters. For 
this reason, and as the observational and physical uncer- 
tainties (slit position, streaming motions) make it diffi- 
cult to determine a physically meaningful rotation veloc- 
ity with an accuracy of more than a few km s~^, we have 



imposed a minimum error on each point of 4 km s~^ (be- 
fore inclination correction). The curves were corrected for 
inclination using the values given in Tables 1 and 2. 

The end result is a smooth representation of the data, 
which is reproducible and as objective as possible, to use 
as input for the mass models. Figure 1 show overlays of 
both the raw hybrid curves and the smooth versions. It 
is easy to verify that no systematic differences in slope or 
shape have been introduced. The error bars in the smooth 
curves are also a good representation of the uncertainties 
in the underlying raw data. 

Table 3 contains the hybrid smooth rotation curves. For 
each galaxy we list the radii in arcsec, as well as in kpc, 
together with the observed rotation velocities and the un- 
certainties in these values. Also included are the rota- 
tion curves for the gas component (already included is the 
factor 1.4 mass scaling for He), the disk component (val- 
ues listed assume T^,{R) — 1.0), and where applicable the 
bulge component (also for T^,{R) = 1.0. 

2.3. Comparing the smooth curves 

As noted above, we have included the five LSB galax- 
ies presented in SMT in our Sample I. As SMT show their 
raw data and derived smooth curves, we can compare both 
sets of smooth rotation curves to investigate possible sys- 
tematics in our respective methods. This is done in Fig. 2. 
It is clear that the correspondence between both velocities 
and errorbars is good and the differences are minor. In 
most cases (F568-1, F568-V1 and F574-1) both sets agree 
at (better than) the la level. Small remaining differences 
are usually caused by a slightly different estimate of the 
velocities in sparsely sampled parts (F563-V2 and F568-3). 
The SMT curve for F563-V2 (Fig. 2, top-left) is shghtly 
higher than our curve in the inner parts, but falls below 
ours in the outer parts. The sparscness of optical data 
points in the outer parts and different interpretations of 
the continuity between Hi and optical data are probably 
the main cause of this difference. For F568-3 (Fig. 2, bot- 
tom) we find a small (< la) systematic difference between 
both smooth curves. This galaxy has been measured in- 
dependently by SMT and by us (Paper I). These raw data 
sets agree in detail, and the difference must therefore be 
due to a slightly different interpretation of the sparse raw 
data. In summary, the smooth curves we present here give 
a good and reproducable representation of the data. 

3. MASS MODELS 

In order to find the signature of the dark halo one needs 
to model the observed rotation curve using a number of 
separate dynamical components, described below. 

3.1. Stellar component 

To model the stellar disk, the i?-band photometry pre- 
sented in de Blok, van der Hulst & Bothun (1995) was 
used. The rotation curve of the disk was computed fol- 
lowing Casertano (1993) and Begeman (1987). The disk 
was assumed to have a vertical sech^ distribution with a 
scale height zq = h/6 (van der Kruit & Searle 1981). The 
rotation curves of the stellar component were resampled 
at the same radii as the smooth curves. We assume T^. 
is constant with radius. While one expects some modest 
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variation in T^, with radius (de Jong 1996), the color gra- 
dients in LSB galaxies tend to be small, so this effect is 
not likely to be significant. 

3.2. Gas disk 

The Hi surface density profiles presented in BMH and 
van der Hulst et al. (1993) were used. They were scaled 
by a factor of 1.4 to take the contribution of helium and 
metals into account. Their rotation curve was derived as- 
suming the gas was distributed in a thin disk. The gas 
rotation curves were resampled at the radii of the smooth 
observed rotation curve. 

3.3. Dark halo 

The dark halo component differs from the previous two 
in that we are interested in parametrizing this component 
assuming some fiducial model. The choice of this model is 
the crux of most of the DM analyses in the literature, and 
many models exist. These can be broadly distinguished 
in two groups: halo models with a core, and halo mod- 
els with a cusp. An example of the first category is the 
pseudo- isothermal halo, an example of the latter the CDM 
NFW halo. 

As one of the goals of this paper is to assess the rele- 
vance of either category to the high-resolution LSB galaxy 
rotation curves, we will present models derived using both 
models. We do realize there are many intermediate models 
described in the literature that probably can fit our data 
equally well. However, our goal here is simply to see where 
the data lead us: is there a preference for models with a 
core or with a cusp? We now describe the details of both 
models. 

3.3.1. Pseudo-isothermal halo 

The spherical pseudo-isothermal halo has a density pro- 
file 

where po is the central density of the halo, and Re the 
core radius of the halo. The corresponding rotation curve 
is given by 



Piso{R) = Po 
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(1) 



l-^arctan(A) 
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V{R) = ^/AirGpoUc 
The asymptotic velocity of the halo, Kx) , is given by 



(2) 



inGpoRc 



(3) 

To characterize this halo only two of the three parameters 
{po,Rc,yoo) are needed, as equation (3) determines the 
value of the third parameter. 

3.3.2. NFW halo 



The NFW mass density distribution takes the form 
Pnfw{R) ^ 



2 (4) 

{R/Rs) (1 + R/Rsf 

where Rs is the characteristic radius of the halo and pi is 
related to the density of the universe at the time of col- 
lapse. This mass distribution gives rise to a halo rotation 
curve 

n 1/2 

ln(l -|- ex) — ex/ (1 + cx^^ 



V{R) = V-2 



200 



_ x[\n{l + c)^cl {I + c)] 



(5) 



where x — -R/i?200- It is characterized by a concentration 
parameter c — i?20o/-Rs and a radius -R200- These are di- 
rectly related to Rg and p^, but are used instead as they 
are a convenient way to parametrize the rotation curve. 
The radius i?200 is the radius where the density contrast 
exceeds 200, roughly the virial radius (Navarro, Frenk & 
White 1996). The characteristic velocity V200 of the halo 
is defined in the same way as i?2oo- These parameters are 
not independent and are set by the cosmology. 



3.4. Mass-to-light ratios and weighting 

One of largest uncertainties in any mass model is the 
value of T^,. Though broad trends in T^, have been mea- 
sured and modelled (e.g. Bottema 1997; Bell & de Jong 
2000a), the precise value for an individual galaxy is not 
well known, and depends on extinction, star formation his- 
tory. Initial Mass Function, etc. Rotation curve fitting is 
a problem with too many free parameters (van Albada & 
Sancisi 1986; Lake & Feinswog 1989) and some assump- 
tions regarding T^, must be made. We therefore present 
disk-halo decompositions using four different assumptions 
for T^ for the galaxies in Sample I. For the galaxies in 
Sample H only the minimum-disk model is presented. 

Minimum disk. This model assumes that the observed 
rotation curve is due entirely to DM. This gives an upper 
limit on how concentrated the dark mass component can 
actually be and is the version of minimum disk preferred 
in the CDM literature. 

Minimum disk + gas. The contribution of the atomic 
gas (Hi and He) is taken into account, but T^ is assumed 
to be zero. This is the classical definition of minimum disk 
as used in the Hi rotation curve literature. 

Constant T^. Here T^ is set equal to a constant 
value based on an Initial Mass Function and a star for- 
mation history appropriate for LSB galaxies. For the 
range in color 0.4 < B ~ V < 0.65 which LSB galax- 
ies normally exhibit (de Blok, van der Hulst & Bothun 
1995) a value Ti,{R) = 1.4 is a good estimate. For ex- 
ample, using the Bruzual & Chariot (1993) model with 
constant star formation rate and Salpeter IMF, we find 
that T*(i?) = 1.4 corresponds to B - V ^ 0.46. The 
PEGASE2 model (Rocca-Volmerangc, priv. comm.) gives 
a value B — V = 0.38, whereas the model by Cole et al. 
(2000) yields B ~V = 0.67. The models by Bell & de 
Jong (2000b) give values around B — V c:^ 0.6. The value 
T^. = 1.4 is thus actually at the "light-weight" end of the 
plausible range, but this was deliberately chosen in order 
to give maximum opportunity for the cuspy NFW models 
to fit the data. We realize that the values derived here 
should not be regarded as definitive: changes in the IMF 
model used or different estimates for internal extinction 
can lead to different values. However, here we attempt 
to derive a conservative estimate for T^, based on the ob- 
served properties of the stellar population. Further (up- 
ward) refinement of the T^. value is thus more likely to 
cause more problems for NFW fits. 

Maximum disk. The rotation curve of the stellar com- 
ponent is scaled to the maximum value allowed by the 
(smooth) rotation curve, but with the restriction that the 
DM density is required to be positive at all radii (thus 
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avoiding a so-called "hollow halo" ) (van Albada & Sancisi 
1986). Because of the different DM distributions that we 
test (core and cusp) this can ocassionally lead to maximum 
disk values that differ slightly for each of the two models. 
A more extensive description is given in §6. 

Each of the rotation curves was fitted using the gipsy 
task ROTMAS. The program determines the best-fitting 
combination of Re and Vco (for the pseudo-isothermal 
halo) or c and V200 (for the NFW halo), using a least 
squares fitting routine. We assigned weights to the data 
points inversely proportional to the square of their uncer- 
tainty. Additional checks were made with other fitting pro- 
grams to check the results of the fits (discussed in Sect. 4.1 
below). We also refitted the smooth rotation curves pre- 
sented in SMT and were able to reproduce the numbers 
given in their Table 2. 

4. RESULTS 

Tables 5 and 6 give the results of the model fitting us- 
ing the NFW halo model. Tables 7 and 8 show the results 
for the pseudo-isothermal halo. Figure 3 presents the re- 
sults of the NFW halo and pseudo-isothermal halo mass 
modelling for each galaxy side-by-side. The two leftmost 
columns show the results for NFW halo fitting. The two 
rightmost columns show the pscudo- isothermal halo fitting 
results. 

The first and third columns show the best-fitting mod- 
els. The rotation curves of the gas are shown as dotted 
lines, those of the stellar disk component as short-dashed 
lines. In two cases (U6614 and F571-8) a significant bulge 
was present which was modelled separately; this is shown 
as the dot-dashed line (see also Section 4.3). The resulting 
halo rotation curve is shown as the long-dashed line. The 
final total model curve is drawn as the full line. In each 
of the model panels we also give the reduced x^ of the fit 
and the chance p that the data and the model could result 
from the same parent distribution. This probability was 
derived using a simple x^ test; it is an indicator for the 
compatibility of the data and the model chosen to describe 
it. Values p > 0.95 indicate that the data and the model 
are a good match. Values p < 0.05 indicate that the model 
is incompatible with the data, and that better models can 
be found. 

We show two best-fitting values. One as found by the 
hnear fit of the gipsy rotmas task (indicated as a "plus"- 
cross (-1-)), and one found by finding the minimum in 
the plotted logarithmic parameter space (indicated by a 
"times" -cross (x)). These two are identical, except when 
extreme parameter values occur (usually during maximum 
disk fits) and numerical precision of the fitting routine 
starts to play a role. This effect is visible in the bottom- 
right corners of the NFW contour plots, where the very 
large values of V200 in combination with the small values 
of c cause increasingly ragged contours. A large differ- 
ence between these two best values therefore indicates that 
the fit should not be regarded as definitive. Indeed, in a 
number of these cases (indicated in the Tables as italic 
numbers) the fitting routine was unable to determine a 
valid solution, and an indicative value had to be chosen 
by hand. This happened mostly with the NFW models. 
One of the reasons for this is that the inner parts of the ro- 
tation curves can be well described by V{R) ~ R, whereas 



the NFW model has the form V ^ R^/^ . To accommodate 
the model, the fit tries to stretch out the NFW curve (re- 
sulting in small c and high V200) in order to make it look 
linear. 

The second and fourth columns of panels in Fig. 3 show 
the Icr (thick contour) and (2,3,4, 5)tT (thin contours) 
probability contours of the halo parameters in logarith- 
mic space. The reason for choosing a logarithmic repre- 
sentation is that the y^ distributions for the NFW halo 
parameters often show extended tails towards very small 
c and/or large V200- For comparison. Fig. 4 shows a repre- 
sentative example of the cr-contours for the minimum disk 
model of F583-4 plotted in linear C-V200 space. 

It is important to realize that the a contours are plotted 
with respect to the minimum x^- That is, existence of a 
narrow distribution only means that the minimum x^ is 
well defined. It does not imply that the fit is good in an 
absolute sense. For that one needs to refer to the value 
of the reduced x^ itself or the probability p that the data 
and model are compatible. There are many cases where 
the NFW model is not a good fit, making it difficult to 
plot absolute likelihood contours. 

Finally, in the NFW contour plots in the second column 
we show the range of c and ^200 values for the currently 
popular ACDM cosmology as derived from numerical mod- 
els (Navarro, Frenk & White 1997) (see §5.3). The cross- 
hatched and single hatched grey areas shows the expected 
1(7 and 2(7 logarithmic scatter in c (where (Jc ~ 0.18) as 
found in numerical models by Bullock et al. (1999). In- 
dependent simulations by Jing (1999) find a much smaller 
logarithmic scatter of cTc = 0.08. The latter do of course 
put much stronger constraints on the NFW results. For 
the sake of clarity, however, and to give the NFW model 
as much chance as possible we adopt the larger estimate 
of the scatter in c of Bullock et al. (1999). 

For the pseudo- isothermal halo we show contours of con- 
stant central density p^. The contours represent from top 
to bottom pq ^ 0.1 (dotted), 1 (dashed), 10, 100, 1000 
(dotted) xl0"3 Mq pc"^. 

4.1. Weighted versus uniform 

To investigate how stable the derived halo parameters 
are with respect to the precise definition of the errorbars, 
we have re-derived the models assigning uniform and equal 
weights to all data points. Though we do not list the latter 
values here, we show in Fig. 5 a comparison between the 
two sets of parameters. It is clear that these agree well, 
showing that the results presented here are robust against 
the precise definition of the errorbars. 

4.2. Mass models: smooth and raw 

Has the procedure used to derive the smooth rotation 
curves affected some of the model results? We established 
in Sect. 2.4 that this procedure introduced no systematic 
differences between the smooth curves and the raw data. 
Here we test this again by checking whether the smooth 
curves give the same fit results as the raw data. 

As the NFW model is more sensitive to changes in the 
inner slope than the pseudo-isothermal model we will use 
the former in our checks. We first fit NFW minimum disk 
models to the smooth and the raw curve of F583-1, as 
a representation of the data from Paper I. Both fits are 
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presented in Fig. 6, where we have imposed a minimum 
error of 4 km s^^ on the raw data to make the errorbars 
consistent with the smooth curve. It is clear that the two 
fits are identical within the errorbars. Similar results are 
obtained using other curves from Paper I. 

As a second test we evaluate the SMT data. We have 
fitted several minimum disk NFW models to each of the 
SMT galaxies. We have fitted the SMT raw data, the 
smooth curve presented in SMT, and our smooth curve 
derived from the raw SMT data. These fits were done in- 
dependently by two of us using independent fit codes on 
the smoothed (dB) and unsmoothed (McG) data. Table 4 
lists the derived parameter values. For comparison, we 
also list the results for our own independent observation 
of F568-3 from Paper I. In Fig. 7 we compare the c values 
derived for each galaxy. 

We see that the galaxies for which our smooth curves 
and the ones presented in SMT agree, also have similar 
model parameters, which agree with those derived from 
the raw data (F574-1 and F568-V1). In the other three 
cases (F563-V2, F568-1, F568-3) the c- values derived from 
our version of the smooth curves agree with those derived 
from the raw data, whereas the SMT c-values are higher. 
It is important to keep in mind that even though the for- 
mal fit values show a large discrepancy, the rotation curves 
themselves only show very subtle differences (Fig. 2). This 
illustrates the importance of having high-accuracy rotation 
curves of a large sample. In the following we only consider 
our smooth versions of the SMT data. 

In summary, we believe that the results from our smooth 
curves are not systematically different from the raw data. 
As stated before, we prefer to use the smooth curves as 
these are more evenly sampled, and prevent the occurence 
of imaginary halo masses which can arise when the oc- 
casional (raw) data point happens to scatter below the 
rotation velocity of the disk alone. 

4.3. Remarks on individual galaxies 

F563-1 For this galaxy independent observations are 
available from de Blok & Bosma (2001). See Paper I for 
a comparison. Note that the observed curve differs sig- 
nificantly from the "beam-smearing corrected" model pre- 
sented in van den Bosch et al. (2000). The model presented 
there shows an almost flat rotation curve over most of the 
radial range, which clearly disagrees with the new data. 
Beam-smearing corrections are not infallible. 

F563-V2 This is our version of one of the SMT curves. 
This curve does significantly worse at fitting NFW than 
a pseudo-isothermal halo. The systematics seen here are 
typical for many of the NFW fits: the inner parts are 
overestimated; the model then underestimates the mid- 
dle parts, and shoots up again in the outer parts. For 
this galaxy no i?-band photometry is available and we 
have used B band photometry from McGaugh & Bothun 
(1994). Assuming B — R = 0.9, which is the typical color 
for an LSB, this yields a value for the constant T^, case of 
Ti,{B) = 1.1. The maximum disk NFW model fits signif- 
icantly better than the other NFW models. It is however 
not compatible with cosmological predictions from the nu- 
merical models. 

F568-1 This is another of the curves presented in SMT. 
The systematics of overestimating the inner part, underes- 



timating the middle, and overestimating the outer veloc- 
ities again are also present here. Again maximum disk is 
the best of the NFW models, which is another way of say- 
ing that the shapes of the (inner) rotation curves are more 
like that expected for the stars (albeit with the wrong T^,). 

F568-3 This is a well determined curve, for which there 
are several consistent independent measurements (see Pa- 
per I). 

F571-8 This is the only edge-on galaxy in Sample I, 
so we are concerned about optical depth and projection 
effects in the optical data. (The Hi data are not used 
for this galaxy.) These effects could cause us to measure 
the rotation velocity at a ring where the optical depth be- 
comes unity. Recently, Matthews & Wood (2000) have 
used radiative transfer models to investigate optical depth 
effects on rotation curves in edge-on LSB galaxies, and 
they conclude that these effect are likely to be small due 
to the low dust content in LSB galaxies. Bosma et al. 
(1992), in a comparison of the optical and Hi curves of 
the edge-on galaxy NlOO, also finds that late-type galaxies 
tend to be transparant, even when seen edge-on. Never- 
theless, we cannot exclude the possibility that the shape 
of the optical curve is affected. If this is so then in this 
case the mass model will change. For the constant T^, and 
maximum disk cases we have added an exponential bulge 
with T^(i?)buigo = 0.5. See BMH for a description of the 
bulge-disk decomposition. 

F574-1 Another SMT curve. While the pseudo- 
isothermal halo fits for this galaxy are good, the NFW 
fits show the by now familiar discrepancy: too steep in 
the inner part, underestimating the middle, and rising too 
quickly in the outer parts. Maximum disk NFW provides 
a good fit, albeit with low c and high T^,. F574-1 was the 
worst case of beam-smearing from the Hi sample, but the 
increase in the initital rate of rise of the rotation curve 
found optically does not really help NFW. The optical 
data imply a 'cusp' slope (p oc r") of a = —0.49 ± 0.26 
(de Blok et al. 2001), still weU short of the NFW value 
a = —1. This is the limit in the minimum disk case; if 
allowance is made for stellar mass, a value even closer to 
a constant density core is required. 

F583-1 A well resolved and well-observed curve that 
shows the NFW over/under/over-fit discrepancy. For all 
assumptions about stellar mass, Xiso ^ Xnfw This 
galaxy strongly prefers a halo with a constant density 
core over one with a cusp, a conclusion which has not 
changed from McGaugh & de Blok (1998). Only a sub- 
stantial change in the shape of the rotation curve would 
alter this conclusions, which would require a large system- 
atic error. Beam smearing can no longer be invoked as the 
cause of such a systematic error now that this object has 
been resolved to sub-kpc scales. 

U5750 This curve was observed both by us and de Blok 
& Bosma (2001), and the two data sets show good agree- 
ment (see Paper I) . This curve is difficult to model with a 
standard NFW profile, but the pseudo-isothermal model 
provides a good fit. The outermost point is taken from 
the Hi curve, and provides an important constraint for the 
NFW model. Without this point the fit produces c ^ 
and an impossibly large V20O) a result of the fitting pro- 
gram trying to make V{R) ~ R^'^ look like V ^ R. 

U6614 This is the only giant LSB in Sample I. The anal- 
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ysis is complicated by the presence of a dominant bulge, 
which we have modelled as an exponential spherical bulge 
with h = 3.0" and HoiR) = 18.4 mag arcsec^^. The disk 
has parameters h = 19" and HoiR) — 21.3 mag arcsec"^. 
The rapid rise and subsequent dip in the rotation velocity 
at small radii clearly suggest the dominance of the bulge 
in this giant LSB galaxy. We have assumed the bulge to 
be maximal at T*(i?) = 3.7. As the bulge accounts for 
most of the rotation velocity in the inner parts, this seems 
to argue against cuspy halos in giant LSB galaxies. 

5. DISCUSSION 
5.1. NEW and pseudo-isothermal: a comparison 

The pseudo-isothermal halos generally provide better 
fits than the NFW halos. In Fig. 8 we compare the re- 
duced x^ values for the four different T^ cases. For the 
minimum disk case we plot Samples I and II, for the other 
cases only Sample I is plotted. It is clear that the large 
majority of the curves presented here is best fitted by a 
pseudo-isothermal halo. This holds true even in the max- 
imum disk case where one might naively expect the dom- 
inance of the optical disk to wipe out any discrepancies 
of a particular halo model (though perhaps not for LSB 
galaxies) . 

Another way of comparing the models is given in Ta- 
ble 9. This lists the number of galaxies in Sample I that 
have good {p > 0.95) or bad (p < 0.05) fits for each of 
the two models. Here again it is clear that the pseudo- 
isothermal model performs much better, for every assump- 
tion of T^. These statements do not depend on the errors. 
If we double (or halve) the size of the error bars, x^ '^iU 
change for both halo cases, but will always remain less for 
the pseudo-isothermal case. To alter this result would re- 
quire systematic changes to the shapes of all the rotation 
curves. 

Fig. 9 shows the residuals of the best-fitting minimum 
disk models versus the observed data. Residuals are plot- 
ted against halo scale size (i?200 for NFW and Re for 
pseudo- isothermal halos), radius in kpc, number of optical 
disk scale lengths and fraction of maximum radius of the 
rotation curve. As described in the previous section, the 
NFW fits that fail do so in a systematic way: the inner 
velocity is overestimated, then the model drops below the 
observed velocities in the middle and in the outer parts it 
once again overpredicts the velocity. The NFW residuals 
are most pronounced when plotted against R/Rmax- The 
majority of the residuals change sign at ^ 0.2i?max and 
^ 0.7i?niax- As the radius i?max does not have any physi- 
cal significance, but is determined by the observations (slit 
angle, presence of Ha, etc.) this indicates that the sys- 
tematics are due to the choice of model, rather than being 
associated with any particular length scale in the galax- 
ies. Similar conclusions are reached when the residuals are 
plotted for the min+gas, constant T* and maximum disk 
cases. 

Though this is not readily apparent in Fig. 3, the resid- 
uals for the pseudo-isothermal halo model also show a sys- 
tematic behaviour, though at a much lower level than the 
NFW model. Here the residuals do not increase towards 
the center, and as the typical size of the residuals is smaller 
than the uncertainty in the individual data points, this just 
shows us that the rotation curve shape is subtly different 



from that of a pure pseudo-isothermal halo. This should 
come as no surprise given the simplifying assumptions of 
e.g. minimum disk that we have made. 

5.2. The pseudo-isothermal halo 

Of the two models investigated, the pseudo-isothermal 
halo best describes the data. Here we briefly explore some 
correlations between the pseudo-isothermal halo model pa- 
rameters and the parameters describing the luminous com- 
ponents of the galaxies. To increase the range of the pa- 
rameters we also consider the samples of Broeils (1992) 
of (mainly) luminous HSB galaxies, and Swaters (1999) 
of late-type dwarf galaxies. From these samples we only 
select bulge- less galaxies brighter than Mb = —16.5 to 
be consistent with the range of luminosities found in our 
sample. 

Fig. 10 presents the results for the three samples. We 
show the min-|-gas case, which the two comparison sam- 
ples refer to as their "minimum disk" . The most obvious 
correlation visible in Fig. 10 is that between Re and po- 
This is a reflection of the fact that these two are correlated 
through the asymptotic velocity Voo (See Eq. 3). Lines of 
constant Voo have a slope of — ^ in the Re — po diagram, 
and the diagram therefore just reflects the limited range in 
Knax in the samples. There is an indication that that Re 
increases towards lower surface brightnesses, and that po 
decreases (as one would expect if LSB galaxies inhabited 
lower density halos). 

The large scatter in these figures sheds little light on 
galaxy formation or the details of pseudo-isothermal ha- 
los. How the observed regularities of galaxy kinematics 
(like the TuUy-Fisher relation) can emerge from this scat- 
ter remains a mystery. 

A further analysis is presented in de Blok et al. (2001) 
where the mass density distributions that give rise to the 
observed rotation curves are presented. They show that 
the minimum disk mass density distributions at small radii 
can be parametrized by a power law p ~ r" where for 
the LSB galaxies a = —0.2 ± 0.2, clearly different from 
— 1.5 < a < —1 as predicted by CDM. These minimum 
disk slopes are upper limits. When stars are properly 
taken into account, assuming some reasonable value for 
T^,, the slopes decrease, and become even more consistent 
with constant density cores. SuccesfuU theories of galaxy 
formation and evolution that attempt to model LSB galax- 
ies should thus be able to produce halos dominated by 
constant-density cores. 

5.3. The NFW halo 

As noted earlier, the c and V200 halo parameters are 
related. Here we compare the derived c and V200 values 
with those predicted by ACDM with the Navarro, Frenk 
& White (1997) prescription for flm = 0.3, f^A = 0.7, and 
h — 0.65 with a COBE normalized power spectrum. The 
values of c depend on the assumed T*, which, as discussed 
before, is uncertain. Minimum disk however gives strong 
upper limits on the values of c: when T^. is increased the 
halo needs to compensate by becoming less concentrated 
(Navarro, Frenk & White 1997). Minimum disk models 
with c values higher than found in simulations can usu- 
ally be reconciled with these observations by increasing 
T^, or introducing a bulge to bring the c-values down, as 
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one can see from the progressive decrease in c values from 
minimum disk to maximum disk. 

Explaining minimum disk models with concentrations 
lower than the simulated values is more difficult, ft in- 
dicates one or more of three problems: failure of the 
model, failure of the assumption of circular motion in de- 
riving rotation curves, or a dramatic (non-cosmological) 
re-distribution of DM. This last option is not really un- 
derstood, and potentially removes any of the predictive 
power that the CDM theory has. We will not discuss it 
here, except to note that the most plausible effect, adi- 
abatic contraction, further concentrates the DM, making 
the problem worse. 

As an aside we note here that the min-|-gas case some- 
times gives slightly higher c-values than the simple min- 
imum disk case. In most cases this is due to a central 
depression in the Hi surface density that gives rises to 
imaginary rotation velocities, which have to be compen- 
sated for by the halo. Also some of the outer rotation 
velocity is explained by the gas rotation curve, yielding a 
halo curve that bends more at small radii. Consequently 
the halo model tends to be slightly more concentrated. 

Fig. II shows the derived c and V200 values and com- 
pares them with the ACDM predictions. For the minimum 
disk case we show both Samples I and II, for the other 
three T* values only Sample I is shown. The data points 
are coded to indicate their significance level p. 

Several points can be made about Fig.II. Firstly, the 
bottom-right plot clearly shows that maximum disk is 
inconsistent with the NFW halos expected to arise in 
ACDM. Secondly, there is some correlation between the 
significance of the fits and their position in the c — V200 
diagram. In the minimum disk case the majority (11 out 
of 14) of the p > 0.95 points are found at V200 ^ 100 km 
s~^. Most (17 out of 19) of the p < 0.95 points are found 
to the right of this line. This division becomes more clear 
in the minimum disk-|-gas and constant T^, plots. As the 
high V200 values tend to occur at lower c, this is likely to 
be the effect of a NFW halo trying to fit a solid-body-like 
curve, by hiding its curvature outside the visible galaxy, 
i.e., by decreasing c and increasing V200- 

Thirdly, the distribution of points does not agree with 
that predicted by the numerical models (Jing 1999; Bul- 
lock et al. 1999). There are more points above, but more 
importantly, below the la lines than expected. This low-c 
tail consists of fits that have a high to reasonable signifi- 
cance p associated with them. We show the distribution 
again in Fig. 12. The two histograms are for minimum disk 
(open histogram. Samples I and II) and constant T^, (filled 
histogram. Sample I). Over-plotted are log- normal distri- 
butions showing the distribution derived from numerical 
simulations. Unfortunately, these simulations do not agree 
on the value of the dispersion. The ACDM model by Bul- 
lock et al. (1999) gives a logarithmic dispersion CTc — 0.18, 
while the distribution for relaxed ACDM halos as found 
by Jing (1999) has a logarithmic dispersion of CTc — 0.08. 
The observed distribution is clearly wider than either the- 
oretical one. By changing the cosmology of the model one 
can change the mean of the distribution (e.g. OCDM has 
a mean logc = 1.25) (Jing 1999), but the width hardly 
changes. Thus one can possibly shift the model to higher 
c to fit the high-c end of the distribution, but it is impossi- 



ble to explain the large observed low-c tail with the kind of 
log-normal distribution one derives from the simulations. 

5.4. Morphology 

Rotation curves have the implicit assumption of circu- 
lar motion. Can non-circular motions affect the rotation 
curves? As the NFW models show the largest residuals 
in the centers of some of the LSB galaxies, it is possible 
that they could be affected by non-circular motions due 
to non-axisymmetric components. We will investigate the 
matter here by comparing the morphology of our galaxies 
with the quality of the fits. 

Table 10 contains a short description of the morphol- 
ogy of the galaxies, where we have focussed on the central 
parts. In the Table "core" refers to a galaxy whose cen- 
tral light distribution can best be described by an axisym- 
metric model, presumably implying negligible non-circular 
motions in the inner part. The word "core" is used very 
loosely here. It does not necessarily indicate the presence 
of a bulge or massive central component, but is just an in- 
dication of the (deprojected) round shape of the isophotes 
in the inner part of the galaxy. "Bar" indicates a cen- 
tral morphology dominated by a bar-like structure, usually 
magellanic, that may indicate the presence of non-circular 
motions. 

The results are summarized in Table 11 (for the min- 
imum disk assumption). The conclusion is that there is 
no clear dependence of residual velocity on morphology. 
There is thus no indication that the failure of NFW to 
fit some galaxy rotation curves can be attributed to the 
presence of bars or non-circular motions. 

6. THE MAXIMUM DISK 

As noted in the Introduction, the inner rotation curves 
of HSB galaxies can usually be well explained by scaling 
up the rotation curve derived from the light distribution. 
This maximum disk procedure results in T^, values that 
are reasonably consistent with those derived from stel- 
lar population synthesis models (Verheijen 1997; Palunas 
& WiUiams 2000; van Albada & Sancisi 1986). Further- 
more, bars seem to demand near-maximal disks in HSB 
galaxies (Debattista & Sellwood 2000; Weiner, Sellwood 
& WiUiams 2000) 

The matter of maximum disk in LSB galaxies was 
first discussed by de Blok & McGaugh (1997) who noted 
that from a stellar population point of view maximum 
disk demanded unreasonably high T^ values. Substantial 
amounts of DM were still needed within the optical disk 
to explain the observed Hi curves. SMT revisited the sub- 
ject and noted that the shghtly steeper slopes they found 
using their Ha curves, enabled them to scale up the disk 
rotation curve by an even larger factor. The maximum 
disk Ti, value is extremely sensitive to the inner slope, 
and only a very small increase is needed to change it by 
a significant factor. In LSB galaxies, the maximum disk 
T^ values are thus much larger than expected on the basis 
of colors, metallicities and star formation histories. The 
higher T^, values as found by SMT worsened this prob- 
lem (a consequence already noted in de Blok & McGaugh 
1997), despite the fact that their maximum disk models 
could slightly better reproduce the observed inner curve. 

Though it seems unlikely that the maximum disk results 
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can be explained using "reasonable" stellar populations, 
given what we know about the star formation history, dust 
content and metallicity of LSB galaxies, the matter is still 
relevant for exploring possible baryonic disk DM scenarios. 
To explain maximum disk in LSB galaxies purely in terms 
of baryons one has to assume a large amount of unseen ma- 
terial in the form of e.g. cold molecular gas, optically thick 
neutral hydrogen, low mass stars or non-standard initial 
mass functions. It should be noted though that many of 
these hypothetical mass components would however vio- 
late constraints imposed by disk stability (Athannasoula, 
Bosma & Papaioannou 1987; Mihos, McGaugh & dc Blok 
1997), near-infra red colors (Bell ct al. 2000; Bell & dc Jong 
2000a), and would possibly introduce a surface brightness 
segregation in the baryonic TF-relation (McGaugh et al. 
2000) 

In Fig. 13 wc compare the maximum disk _B-band T^ ra- 
tios^ for Sample I, with those derived by Broeils (1992) for 
a sample of mostly HSB galaxies and by Swaters (1999) for 
a sample of dwarfs. We again show only bulge-less galax- 
ies brighter than Mb = —16.5. Also indicated are T* 
values from Bell & de Jong (2000b), who tabulate stellar 
mass-to- light ratios for various star formation histories and 
population synthesis models as a function of color. Here 
we show representative values (assuming a simple Salpeter 
IMF) spanning the color range exhibited by late-type HSB 
galaxies, gas-rich dwarfs and LSB galaxies. 

The T^, values for HSB galaxies agree to within a factor 
of 3 and can be considered to be (close to) maximum disk. 
The values found for LSB galaxies and dwarfs are less eas- 
ily reconciled with the model values. Observationally val- 
ues up to Ti,{B) = 15 are found, while the typical model 
value (again for a simple Salpeter IMF) is T^,{B) ~ 0.9 
ioi B — R = 0.8 (the average color for a dwarf/LSB). This 
discrepancy cannot be explained with extinction or pop- 
ulation effects. Extinction in dwarfs and LSB galaxies is 
less than in HSB galaxies (TuUy & Verheijen 1998), and 
a factor of ^ 17 (3.0 mag) extinction is hard to reconcile 
with the known properties of LSB galaxies. Line-of-sight 
extinctions observed towards Hii regions in LSB galaxies 
are never as large (McGaugh 1994; de Blok & van der 
Hulst 1998). 

Apart from changing the IMF in an ad hoc way it is hard 
to see how such high T* values can be reached given the 
constraints imposed by what we know about the star for- 
mation history (low SFR in past and at present) and the 
blue (optical and NIR) colors of LSB galaxies (de Blok, 
van der Hulst & Bothun 1995; McGaugh & Bothun 1994; 
van den Hoek et al. 2000; BeU et al. 2000; BeU & de Jong 
2000a; de Jong 1996). It is likely that the maximum disk 
values as found in LSB galaxies are not representative of 
the evolutionary stage of these galaxies. While the maxi- 
mum disk prescription now has somewhat greater success 
in predicting the inner shape of the rotation curves of LSB 
galaxies, it requires stellar mass-to-light ratios which are 
too large for the stellar populations in these galaxies. The 
mass discrepancies are still large; all this does is move the 
DM from halo to disk. 



6.1. The maximum surface density of a disk 

Just as the minimum disk assumption gives us an upper 
limit on the amount of DM implied by rotation curves, 
the maximum disk hypothesis gives us an upper limit on 
the amount of mass that could potentially be hidden in a 
disk. It is therefore still useful to ask ourselves what these 
maximum disk upper limits imply for the stellar disks. 

Maximum disk means maximum surface density (lumi- 
nous surface density times T^), and therefore gives an ab- 
solute upper limit on the mass surface density in stellar 
disks (for mass components that are distributed like the 
stars). Figure 14 summarizes the maximum disk results 
for the Sample I LSB galaxies as well as the Broeils (1992) 
and Swaters (1999) HSB and dwarf samples. We plot 
the maximum disk T^, values, as well as the luminosity, 
rotation velocity, surface brightness, and maximum disk 
surface density. The data are divided into three surface 
brightness bins. As already shown in de Blok & McGaugh 
(1997), at fixed Vmax LSB galaxies have higher maximum 
Ty, values than HSB galaxies. 

Figure 14 also shows the maximum surface density a. 
As the decrease in surface brightness is faster than the 
increase in maximum disk T^,, towards low surface bright- 
nesses the maximum surface density tr in a disk decreases 
with surface brightness. The no — a panel suggests that 
there is a well-defined upper limit to the maximum surface 
density disks can attain. Even under maximum disk, LSB 
galaxy disks have on average lower surface densities than 
HSB galaxy disks, again putting limits on the amount of 
baryonic mass one can hide in these disk. 

Figure 3 shows that even in the maximum disk case most 
LSB galaxies have Vmax(disk) < Vmax (observed). There- 
fore, even in the maximum disk case a moderate amount 
of DM is still required in the optical disk. It is therefore 
hard to explain the TF relation for LSB galaxies in the 
context of maximum disk: the stellar disk then needs to 
provide the luminosity, and the necessary rotation veloc- 
ity. LSB galaxies would deviate systematically from the 
HSB TF relation, which is evidently not the case (Zwaan 
et al. 1995). 

This is illustrated in the inset panel in Fig. 14 (lower- 
right). Using the arguments in Zwaan et al. (1995) we find 
that So(T)^ needs to be constant for galaxies to obey a 
TF relation independent of surface brightness. If all galax- 
ies were truly maximum disk (in the sense that 14nax(disk) 
— Vinax (observed) , one could replace this by the require- 
ment that So(T^;)^ax needs to be constant. The lower- 
right panel shows that this is not the case: at fixed Vmax 
there is a substantial scatter which would translate in ^^ 5 
mag scatter in TF. Clearly the observed scatter is much 
smaller, and this shows the clear need for an additional 
mass component to make TF work. In other words, max- 
imum disk for all galaxies and TF are incompatible. 

7. CONCLUSIONS 

The most important conclusion from this work is 
that the large majority of the high-resolution rotation 
curves presented here prefer the pseudo-isothermal core- 
dominated halo model. For a small number of galaxies 

^ These were derived by converting our i?-band values using the B — R color. We converted our data to B-band, rather than converting the 
HSB data to i?-band, as the colors of LSB galaxies are better determined than those of the Broeils HSB galaxies. The color gradients in LSB 
galaxies are small, so systematic effects are negligible. 
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neither the pseudo-isothermal nor the NFW models are an 
adequate description of the data. This should not come as 
a surprise as the true DM distribution is likely to be more 
complex than the models presented here. Nevertheless, the 
general trend is that for almost all galaxies discussed here 
the relative quality of the fits using the pseudo-isothermal 
model is better than those for the NFW model. 

For a small numer of galaxies the NFW model provides 
a good fit, but generally the concentrations derived from 
the observed rotation curves are lower than predicted by 
the simulations. This is hard to fix: the most likely ef- 
fect which may alter the initial cosmological NFW halo 
is adiabatic contraction, but this has the effect of making 
the final (observed) halo more concentrated, so one would 
have to start off with (cosmologically relevant) halos that 
are even less concentrated. 

It is worrying that for one or two extreme cases the 
difference between "CDM does work" or "CDM doesn't 
work" depends on subtle differences in data, data han- 
dling, or analysis. Figure 7 illustrates that opposite claims 
can sometimes be made from the same data. Hence we re- 
iterate the need for the highest quality data of a large 
sample, in order to minimize these effects. 

We refer to de Blok et al. (2001) where it is shown 
that all data presented here are consistent with a core- 
dominated model; the good NFW fits that are found for 
a number of LSB galaxies can be attributed to resolution 
effects. 

We summarize our results as follows. 

• Pseudo-Isothermal halos are a better description of 
the data than NFW halos. 

• The number of galaxies that cannot be fit 
with NFW halos is significantly larger than the 
number of galaxies that cannot be fitted with the 



pseudo- isothermal model. 

• The quality of the fit is not obviously related to 
morphology, luminosity, or surface brightness. 

• A larger number of low-c NFW halos is found than 
one would expect based on the distribution derived 
from CDM simulations. 

• If one were to construct models that would have 
the correct C-V200 values as predicted by cosmology, 
then the resulting T^, values would be too low to 
be consistent with stellar population numbers. 
The shape of the curves would still be wrong. 

• The maximum disk prescription works to predict 
the inner rotation curve shape to some extent, but 
gives mass-to-light ratios which are too high to 
be accounted for by stellar population synthesis 
models. 

• Applying the maximum disk values yields absolute 
upper limits on the disk mass surface density that 
is strongly correlated with surface brightness. 
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Fig. 1. — Comparison of the raw hybrid rotation curves (black dots) with the smooth curves (full lines). The derived uncertainties in the 
smooth curves arc indicated by the dotted lines. 
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Fig. 2. — Comparison of our analysis of the STM data with their resanipled rotation curves. The top four panels and the bottom left panel 
show the raw data from STM (grey small dots), their resampled and smoothed rotation curves (open circles), and our local regression fits to 
the same data (filled large circles). The raw data have been offset by +0.2" to avoid overlap with the binned data. The bottom right panel 
shows the raw data for F568-3 taken from Paper I, along with the STM model (identical to the model shown in the bottom-left panel) and 
our resampled rotation curve based on the data from Paper I. Small differences between the various curves are discussed in the text. 



Fig. 3. — Mass models assuming NFW halo (left) and pseudo-isothermal halo (right). Within each panel the left column shows the best 
fitting model, the right column shows the probability distribution of the halo parameters. For a full description see text (Sect. 4). Figure 3 
not available due to size limitations. All figures available at http://MVWMV.atnf.csiro.au/~edeblok/papers/deblok.ps.gz 
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Fig. 4. — The error contours for the minimum disk NFW fit of F583-4 drawn in linear c — V200 space. Contour values are as in Fig. 3. 
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Fig. 5. — Comparison of the NFW halo parameters c and V200 and the pseudo-isothermal halo parameters Re and po, as derived using a 
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Fig. 6. — Comparison of NFW minimum disk fits to the raw rotation curve of F583-1 (left panel) and the smooth curve (middle panel). 
The right panel compares the two fits, which arc virtually identical. The fit parameters shown in the panels agree within their errors. 
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While the formal fits differ significantly in these cases, the differences in the curves being fit are subtle (see Fig. 2). This illustrates the 
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Fig. 8. — Comparison of the reduced x2 values using NFW and ISO halos, using the four assumptions for T* as described in the text. Note 
that the axes have logarithmic scales. The dotted lines are lines of equality. 
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Fig. 10. — Correlations involving pseudo-isothermal halo parameters assuming minimum disk+gas. Shown are maximum rotation velocity 
Vmax [km s^-*^], surface brightness fio[B-mag arcsec"^], central halo density po [10~^ Mq pc~^], halo core radius Rq [kpc], and the ratio 
of core radius and optical disk scale length Rc/h. Included are the bulge-less HSB galaxies brighter than Mb = —16.5 from Broeils (1992) 
and the high-quality (g = or 1) curves of bright Mg < —16.5 dwarfs from Swaters (1999). Black circles: fH}{B) < 21.9 mag arcsec^^; open 
squares: 21.9 < ^o(^) < 23.2; stars: tJ-o{B) > 23.2 mag arcsec"^. See text for details. 
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Fig. 11. — The NFW halo concentration parameter c plotted against the halo rotation velocity V200 for the four different T*^ cases discussed 
in this paper. Black dots represent good fits (p > 0.95), open dots average quality fits (0.05 < p < 0.95). Crosses represent bad fits (p < 0.05). 
Good fits are primarily found at V200 < 100 km s~^. Maximum disk is clearly inconsistent with NFW. The line labeled "ACDM" shows the 
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Fig. 12. — Distribution of the c-parameter for the minimum disk case (open histogram) and the constant T* case (solid histogram). Over- 
plotted is the theoretical log-normal distribution for a ACDM cosmology, derived from independent numerical simulations by Jing (1999) and 
Bullock et al. (1999). The former finds a log-normal distribution with a logarithmic dispersion Cc = 0.08. The latter finds a wider log-normal 
distribution with ac = 0.18. The observed low-c tail is not consistent with either theoretical distribution. The theoretical distributions have 
been arbitrarily normalized to coincide with the maximum of the observed distribution. 
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[Mq/Lq g], maximum observed rotation velocity Vmax [km s~^] and luminosity Mg [mag]. Included are again the bulge-less bright HSB 
galaxies from Broeils (1992) and the high-quality curves of bright dwarfs from Swaters (1999), as described in the previous figure caption. 
Black circles: ^o(B) < 21.9 mag arcsec"^; open circles: 21.9 < fJ.o{B) < 23.2; stars: l^oiB) > 23.2 mag arcsec"^. The inset panel in the 
lower-right corner shows the product Eo(T,)^, where Sq is the central surface brightness expressed in Lq pc~^. This product should be 
constant with small scatter for a maximum disk interpretation of the TF relation. The TF relation is shown in the top-left corner. It has a 
slope of —8.4. The dotted lines represent the Icr and 2cr scatter where a = 0.81 mag. This scatter is reduced to 0.51 mag when the 4 most 
outlying points are omitted (all galaxies with low inclinations). 
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Table 1 
Sample I: galaxies with photometry 



Name 


D 


Mo(B) 


h 


Mabs(S) 


-'i-niax 


*^max 


Vhcl 


i 




(Mpc) 


(mag/D") 


(kpc) 


(mag) 


(kpc) 


(km s-i) 


(km s-i) 


n 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


F563-1 


45 


23.6 


2.8 


-17.3 


17.7 


112 


3502 


25 


F563-V2 


61 


22.1 


2.1 


-18.2 


9.2 


118 


4312 


29 


F568-1 


85 


23.8 


5.3 


-18.1 


14.9 


142 


6524 


26 


F568-3 


77 


23.1 


4.0 


-18.3 


16.5 


105 


5913 


40 


F568-V1 


80 


23.3 


3.2 


-17.9 


19.0 


118 


5768 


40 


F571-8 


48 


23.9^ 


5.2 


-17.6^ 


15.6 


144 


3768 


90 


F574-1 


96 


23.3^ 


4.3 


-18.4^ 


15.4 


100 


6889 


65 


F579-V1 


85 


22.8^ 


5.1 


-18.8^ 


17.3 


114 


6305 


26 


F583-1 


32 


24.1 


1.6 


-16.5 


14.6 


87 


2264 


63 


F583-4 


49 


23.8^ 


2.7 


-16.9^ 


10.0 


70 


3617 


55 


U5750 


56 


23.5^ 


5.6 


-18.7^ 


21.8 


79 


4177 


64 


U6614 


85 


23.4 


8.1 


-20.3 


62 


204 


6371 


36 



Note. — (2) Distance computed assuming Hubble flow after correction for galactic rota- 
tion and Virgocentric flow. (6) Maximum radius of rotation curve. (7) Maximum velocity 
in rotation curve. (8) Heliocentric systemic velocity. Photometric and distance data from 
de Blok & McGaugh (1997). i?max, ^ax and l^yg derived from new optical curves. 

^Converted from i?-band assuming B — R = 0.9. 



Table 2 
Sample II: galaxies without photometry 



Name 


D 


Vhcl 


Mabs(S) 


^max 


»^max 


i 




(Mpc) 


(km s-i) 


(mag) 


(kpc) 


(km s-i) 


11 
(7) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


F730-V1 


144 


10714 




11.9 


145 


51) 


U4115 


3.2 


343 


-12.4 


1.0 


40 


74 


U11454 


91 


6628 


-18.6^ 


11.9 


152 


64 


U11557 


22 


1390 


-20.0 


6.2 


95 


36 


U11583 


5 


128 


-14.0^ 


1.5 


36 


83 


U11616 


73 


5244 


-20.3^ 


9.6 


143 


60 


U11648 


48 


3350 


-21.0^ 


12.7 


145 


90 


U11748 


73 


5265 


-22.9^ 


21.0 


242 


78 


U11819 


60 


4261 


-20.3^ 


11.7 


153 


66 


E0140040 


212 


16064 


-21.6 


29.2 


263 


35 


E0840411 


80 


6200 


-18.1 


8.9 


61 


90 


E1200211 


15 


1314 


-15.6 


3.5 


25 


70 


E1870510 


18 


1410 


-16.5 


3.0 


40 


58 


E2060140 


60 


4704 


-19.2 


11.6 


118 


39 


E3020120 


69 


5311 


-19.1 


11.0 


86 


55 


E3050090 


11 


1019 


-17.3 


4.8 


54 


53 


E4250180 


86 


6637 


-20.5 


14.4 


145 


33 


E4880049 


22 


1800 


-16.8 


6.0 


97 


63 



Note. — Columns (2) and (3): distance D was calculated from 
Vhcl after correcting for galactic rotation and assuming pure Hubble 
flow with i^o = 75 km s^^ Mpc^^. Column (4): Absolute magnitude 
computed using apparent magnitudes from ESO-LV and RC3, and are 
corrected for foreground Galactic extinction. 

'^The apparent magnitude is Zwicky magnitude 17, and therefore 
very uncertain. 
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Table 3 
Modeled Rotation Curves 



Galaxy 


R 


R 


Vgas 


Vd 


isk 


H^ige"'" 


Vobs 


o-y 




(arcsec) 


(kpc) 


(km/s) 


(km/s) 


(km/s) 


(km/s) 


(km/s) 


F583-1 






















0.3 


0.1 


-0.1 




0.4 







1.1 


11.1 




2.8 


0.4 


-0.9 




4.0 







10.0 


7.0 




5.0 


0.7 


-1.5 




6.7 







17.4 


9.6 




6.9 


1.0 


-2.2 




8.5 







23.5 


11.2 




9.0 


1.4 


-2.8 




10.1 







31.0 


5.2 



'^Only given when known (Sample I). Set to zero if unknown. 
•^Assumes Mgas = IAMri- 
'^For M/L = 1.0 in the i?-band. 

Note. — The complete version of this table is in the electronic edition of 
the Journal. The printed edition contains only a sample. These data are also 
available in electronic format from http://www.atnf.csiro.au/^cdcblok/data and 
http://www. astro. umd.cdu/^ssm/data. 









Table 4 












Comparison i 


OF FITTING parameters 








Obs 


Curve 


NFW halo, minimum disk 


Galaxy 


c 


Ac 


V200 


AV200 


v2 


F563-V2 
F563-V2 
F563-V2 


SM'l' 
SMT 
SMT 


SMI' 

dBMR 

data 


16.2 
7.5 
5.9 


3.4 
3.7 
2.2 


84.5 
153.1 
192.3 


10.4 
70.5 
76.4 


2.516 

1.391 

3.42 


F568-1 
F568-1 
F568-1 


SMT 
SMT 
SMT 


SMT 

dBMR 

data 


13.4 
6.4 
8.3 


1.1 
2.5 
1.1 


112.1 
194.6 
154.5 


6.3 
71.9 
18.5 


0.265 

0.804 

3.49 


F568-3 
F568-3 
F568-3 


SMT 
SMT 
SMT 


SMT 

dBMR 

data 


5.1 
1.17 
1.71 


2.9 
0.5 


160.3 
591.0 
400.4 


88.0 
94.8 


2.147 

3.551 

13.9 


F568-3 
F568-3 


Paper I 
Paper I 


dBMR 

data"" 


3.2 
4.6 


3.7 
0.5 


214.6 
168.4 


233.9 

17.7 


2.239 

8.01 


F568-V1 
F568-V1 
F568-V1 


SMT 
SMT 
SMT 


SMT 

dBMR 

data 


14.2 
14.6 
15.8 


0.7 
1.2 
1.1 


91.5 
92.1 

85.7 


2.3 

4.9 
3.8 


0.239 

0.197 

12.7 


F574-1 
F574-1 
F574-1 


SMT 
SMT 
SMT 


SMT 

dBMR 

data 


9.4 
8.3 

8.2 


0.7 
1.3 
0.4 


91.2 
98.3 
99.3 


4.3 

10.4 

3.4 


0.421 

1.595 

3.84 



Note. — Italics indicate estimates, not actual fits. V200 is in km s""'^. 
The column labeled "Obs" gives the source of the raw data. The column 
labeled "Curve" gives the source for the derived rotation curve: "SMT" 
indicates smooth rotation curve from SMT; "dBMR" indicates smooth 
rotation curve from this paper; "data" indicates a fit to the raw data. 

''Uncertain, depends on initial estimates of fit. 
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Table 5 
Fitting parameters NFW halo, Sample I 









minimum disk 








minimum 


disk + 


gas 






Galaxy 


c 


Ac 


V200 


Ay 


v2 


P 


c 


Ac 


V200 


Ay 


Ared 


P 




F563-1 


10.7 


1.2 


93.1 


4.3 


0.092 


0.999 


11.3 


1.3 


87.5 


3.9 


0.089 


0.999 




F568-3 


3.2 


3.7 


214.6 


233.9 


2.239 


0.017 


4.2 


3.3 


161.3 


118.9 


2.386 


0.011 




F571-8 


7.8 


1.1 


163.8 


20.2 


1.501 


0.123 


7.8 


1.1 


163.3 


20.1 


1.477 


0.132 




F579-V1 


20.9 


1.5 


78.4 


2.6 


0.211 


0.998 


22.1 


1.6 


75.1 


2.5 


0.217 


0.998 




F583-1 


5.1 


1.0 


106.6 


17.0 


0.740 


0.746 


6.2 


1.1 


86.6 


12.4 


0.827 


0.648 




F583-4 


5.7 


1.4 


89.5 


19.0 


0.322 


0.944 


6.1 


1.4 


82.2 


15.5 


0.321 


0.945 




U5750 


2.6 


1.5 


123.1 


58.8 


1.243 


0.262 


2.9 


1.6 


105.8 


47.7 


1.203 


0.288 




U6614 


10.3 


2.0 


169.8 


17.7 


4.626 


0.000 


11.0 


2.2 


163.9 


17.1 


4.712 


0.000 




SMT data, 


our analysis 
























F563-V2 


7.5 


3.7 


153.1 


70.5 


1.391 


0.195 


8.3 


3.8 


133.1 


54.1 


1.484 


0.157 




F568-1 


6.4 


2.4 


194.6 


71.9 


0.804 


0.625 


7.5 


2.4 


160.1 


46.7 


0.869 


0.562 




F568-3 


1.2 




591.1 




3.551 


0.000 


1.2 




552.6 




3.573 


0.000 




F568-V1 


14.6 


1.2 


92.0 


4.9 


0.197 


0.999 


15.7 


1.5 


87.2 


5.0 


0.242 


0.997 




F574-1 


8.3 


1.3 


98.3 


10.4 


1.595 


0.085 


9.0 


1.5 


91.1 


9.4 


1.806 


0.041 








constant T*(_R) = 


1.4 








maximum disk 








Galaxy 


c 


Ac 


V200 


AV 


2 


P 


c 


Ac 


V200 


Ay 


v2 
■^red 


P 


j-R 


F563-1 


9.9 


1.2 


88.8 


4.6 


0.089 


0.999 


4.0 


1.2 


110.0 


20.1 


0.098 


0.999 


6.9 


F568-3 


2.3 


4.9 


218.6 


410.9 


2.127 


0.024 


0.4 


19.5 


595.0 


00 


2.015 


0.024 


2.2 


F571-8 


1.6 


5.7 


591.4 




3.776 


0.012 


1.0 




500.0 




7.060 


0.000 


4.2 


F579-V1 


23.2 


1.9 


67.9 


2.4 


0.215 


0.998 


43.4 


14.6 


31.6 


3.9 


0.671 


0.781 


7.9 


F583-1 


5.4 


1.1 


90.6 


14.7 


0.767 


0.716 


2.2 


1.5 


145.5 


78.7 


0.680 


0.804 


6.5 


F583-4 


5.9 


1.4 


76.0 


14.5 


0.271 


0.965 


11.4 


6.1 


12.9 


3.3 


0.196 


0.986 


9.6 


U5750 


1.9 


1.7 


116.9 


80.4 


1.105 


0.354 


1.9 


1.7 


116.9 


80.4 


1.105 


0.354 


1.4 


U6614 


1.7 


1.8 


303.0 


206.5 


4.005 


0.001 


0.4 




145.6 




7.828 


0.737 


7.7 


SMT data, 


our analysis 
























F563-V2 


6.5 


3.8 


136.0 


73.3 


1.047 


0.397 


1.0 




350.0 




0.449 


0.878 


4.1 


F568-1 


6.7 


2.5 


163.8 


56.1 


0.803 


0.626 


0.6 




669.2 




0.636 


0.898 


9.0 


F568-3 


1.0 




519.4 




3.595 


0.000 


1.0 




467.8 




3.624 


0.000 


1.8 


F568-V1 


14.6 


1.4 


88.8 


5.3 


0.228 


0.998 


3.3 


1.5 


125.8 


45.7 


0.222 


0.998 


14.0 


F574-1 


8.2 


1.4 


87.1 


9.5 


1.391 


0.162 


1.5 


1.6 


69.7 


43.9 


0.204 


0.998 


8.1 



Note. — Italics indicate estimates, not actual fits. V2V10 is in km s ^ . 
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Table 6 
Fitting parameters NFW halo, Sample II 









minimum disli 






Galaxy 


c 


Ac 


V200 


AV 


Xred 


P 


F730-V1 


11.8 


1.9 


131.4 


16.2 


0.995 


0.426 


U4115 


5.0 




133.4 




0.777 


0.591 


U11454 


10.4 


2.0 


152.6 


23.3 


3.334 


0.000 


U11557 


1.0 




425.1 




1.367 


0.093 


U11583 


5.0 




93.3 




0.676 


0.641 


U11616 


12.7 


1.8 


124.4 


14.3 


1.254 


0.244 


U11648 


8.0 


0.7 


146.2 


10.9 


0.964 


0.498 


U11748 


52.5 


4.5 


125.2 


3.8 


3.325 


0.000 


U11819 


6.4 


1.9 


252.9 


77.6 


1.348 


0.177 


E0140040 


16.8 


0.8 


203.3 


5.8 


0.152 


0.989 


E0840411 


1.0 




181.0 




1.608 


0.077 


E1200211 


6.4 


2.1 


27.5 


6.9 


0.246 


0.996 


E1870510 


3.8 


1.5 


93.6 


38.6 


0.059 


1.000 


E2060140 


15.2 


1.3 


92.7 


4.4 


0.425 


0.962 


E3020120 


6.6 


1.5 


98.5 


18.4 


0.333 


0.965 


E3050090 


1.0 




323.6 




0.208 


0.999 


E4250180 


3.1 


0.9 


301.6 


77.8 


0.014 


1.000 


E488-049 


4.9 


1.9 


209.6 


90.1 


0.170 


0.999 


Note. — Italics indicate 


estimates, not actual fits. V200 


is in km s~^ 
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Table 7 
Fitting parameters ISO halo, Sample I 









minimum disk 








minimum 


disk + g 


;as 






Galaxy 


Re 


Ai? 


Po 


Ap 


Xled 


P 


Re 


Ai? 


Po 


Ap 


Xred 


P 




^'563-1 


1.72 


0.23 


91.9 


21.6 


0.085 


1.000 


1.55 


0.22 


102.0 


25.2 


0.078 


1.000 




F568-3 


2.92 


0.36 


36.6 


5.4 


0.522 


0.860 


2.71 


0.39 


38.3 


6.8 


0.676 


0.731 




F571-8 


2 12 


0.19 


106.9 


14.0 


1.525 


0.114 


2.12 


0.19 


106.3 


14.0 


1.512 


0.118 




F579-V1 


0.67 


0.02 


574.8 


37.3 


0.026 


1.000 


0.63 


0.03 


630.8 


50.2 


0.037 


1.000 




F583-1 


2.44 


0.06 


33.0 


1.1 


0.037 


1.000 


2.08 


0.11 


37.7 


2.5 


0.117 


1.000 




F583-4 


1.10 


0.13 


85.5 


15.8 


0.329 


0.941 


1.06 


0.11 


88.2 


15.1 


0.267 


0.967 




U5750 


4.25 


0.39 


10.6 


1.0 


0.154 


0.998 


3.96 


0.49 


10.4 


1.3 


0.214 


0.993 




U6614 


1.86 


0.49 


218.4 


102.5 


1.942 


0.021 


1.73 


0.44 


244.7 


111.6 


1.777 


0.040 




SMT data, 


our analysis 
























F563-V2 


1.69 


0.17 


131.2 


19.4 


0.283 


0.972 


1.58 


0.20 


135.3 


24.7 


0.393 


0.925 




F568-1 


2.22 


0.10 


97.2 


6.2 


0.066 


1.000 


2.03 


0.12 


104.6 


8.9 


0.106 


1.000 




F568-3 


3.93 


0.75 


30.2 


5.6 


1.256 


0.245 


3.75 


0.78 


30.8 


6.3 


1.383 


0.173 




F568-V1 


1.45 


0.11 


153.0 


18.2 


0.110 


1.000 


1.33 


0.14 


170.0 


27.8 


0.172 


1.000 




F574-1 


1.83 


0.06 


71.5 


3.7 


0.100 


1.000 


1.70 


0.09 


76.9 


6.3 


0.232 


0.997 








constant T 


^{R) = 


1.4 








maximum disk 








Galaxy 


Re 


Ai? 


Pq 


Ap 


Xred 


P 


Re 


Ai? 


Po 


Ap 


2 
X.red 


P 


JR 


F563-1 


1.72 


0.26 


79.0 


21.3 


0.083 


1.000 


4.09 


1.01 


13.2 


4.9 


0.124 


0.998 


6.9 


F568-3 


3.07 


0.63 


25.7 


6.3 


0.793 


0.623 


3.36 


0.88 


19.7 


5.9 


0.89 


0.536 


2.2 


F571-8 


4.19 


0.28 


34.4 


2.7 


0.405 


0.954 


9.97 


3.43 


9.8 


2.3 


2.639 


0.002 


4.2 


F579-V1 


0.55 


0.04 


694.4 


84.0 


0.06 


1.000 


0.17 


0.16 


1970 


3234 


1.032 


0.415 


7.9 


F583-1 


2.26 


0.11 


31.5 


2.0 


0.103 


1.000 


3.41 


0.24 


14.3 


1.2 


0.121 


1.000 


6.5 


F583-4 


1.02 


0.12 


80.6 


15.2 


0.249 


0.972 


0.23 


0.15 


103.4 


114.1 


0.188 


0.988 


9.6 


U5750 


4.67 


0.74 


7.1 


1.1 


0.262 


0.984 


4.67 


0.74 


7.1 


1.1 


0.262 


0.984 


1.4 


U6614 


12.18 


2.87 


6.3 


1.9 


1.938 


0.022 


112.0 


506.7 


0.4 


0.4 


4.762 


0.000 


7.4 


SMT data, 


our analysis 
























F563-V2 


1.70 


0.24 


96.6 


19.8 


0.298 


0.967 


2.32 


0.63 


30.8 


11.5 


0.191 


0.992 


4.1 


F568-1 


2.11 


0.15 


90.7 


9.0 


0.120 


1.000 


3.08 


0.73 


27.6 


8.5 


0.270 


0.988 


9.0 


F568-3 


4.35 


1.31 


21.5 


6.0 


1.709 


0.065 


4.54 


1.52 


19.4 


5.9 


1.815 


0.046 


1.8 


F568-V1 


1.41 


0.15 


146.1 


23.0 


0.153 


1.000 


3.80 


0.58 


14.8 


2.6 


0.120 


1.000 


14.0 


F574-1 


1.74 


0.09 


63.1 


5.4 


0.182 


0.999 


3.30 


0.83 


4.6 


1.6 


0.108 


1.000 


8.1 



Note. — Re is in kpc. po is expressed in units of 10 "^ Mq pc 
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Table 8 
Fitting parameters ISO halo, Sample II 









inininium disk 




Galaxy 


Re 


AR 


Po 


Ap 


Xled 


P 


^730- VI 


1.46 


0.06 


215.8 


14.2 


0.097 


0.997 


U4115 


0.94 


0.03 


148.2 


2.4 


0.004 


1.000 


U11454 


1.95 


0.10 


146.9 


11.8 


0.423 


0.927 


U11557 


5.48 


0.55 


15.2 


0.9 


0.052 


1.000 


U11583 


0.63 


0.08 


117.8 


16.5 


0.103 


0.999 


U11616 


1.45 


0.05 


208.9 


11.6 


0.140 


1.000 


U11648 


1.95 


0.25 


104.9 


20.7 


3.792 


0.000 


U11748 


0.36 


0.15 


8540 


6661 


5.402 


0.000 


U11819 


2.93 


0.14 


88.2 


5.2 


0.303 


0.991 


E0140040 


2.55 


0.18 


249.5 


27.0 


0.176 


0.982 


E0840411 


6.41 


0.56 


5.2 


0.3 


0.067 


0.999 


E1200211 


0.57 


0.08 


45.5 


9.2 


0.082 


1.000 


E1870510 


0.97 


0.05 


53.5 


3.2 


0.028 


1.000 


E2060140 


1.17 


0.05 


231.1 


16.9 


0.106 


1.000 


E3020120 


1.90 


0.09 


53.6 


3.4 


0.035 


1.000 


E3050090 


2.09 


0.14 


27.3 


1.8 


0.048 


1.000 


E4250180 


4.41 


0.75 


30.0 


6.6 


0.088 


0.997 


E488-049 


1.63 


0.04 


101.1 


3.1 


0.016 


1.000 



Note. — Re is in kpc. 
10^3 M^ pc-3 



Pq is expressed in units of 



Table 9 
Comparison probabilities NFW vs. ISO Sample I 



Number of galaxies (out of 13) 
ISO halo NFW halo 



T. 



p > 0.95 p < 0.05 p > 0.95 p < 0.05 



~8~ 

8 

10 



"T" 
4 
4 
6 



"3~ 
4 
4 
3 



nunimum 
min+gas 
eonstant 
maximum 



Note. — pis the probability that the model is compatible 
with the data. 



LSB rotation curves 



31 



Table 10 
Morphology and probability NFW halos 



Galaxy 


prob.^ 


core^ 


bar"^ 


morphology 


F563-1 


+ 


— 


+ 


mag. irr. 


F563-V2 





- 


+ 


mag. bar 


F568-1 





+ 




spiral 


F568-3 


- 


- 


+ 


spiral with mag. bar 


F568-V1 


+ 


+ 


- 


spiral 


F571-8 





+ 


? 


edge-on 
disk 


F574-1 


- 


+ 


- 


F579-V1 


+ 


+ 


- 


core, flocc. arms 


F583-1 





+ 


~ 


mag. irr. 


F583-4 


+ 


+? 


+? 


fuzzy 


F730-V1 





+ 


- 


spiral 


U4115 





+? 


- 


fuzzy 


U5750 





- 


+ 


mag. bar 


U6614 


- 


+ 


- 


faint, with bulge 


U11454 


- 


+ 


- 


fuzzy spiral, small core 


U11557 





+ 


- 


fuzzy spiral, small core 


U11583 





- 


+ 


faint mag. bar 


U11616 





+ 


- 


fuzzy, irr 


U11648 





- 


+ 


irr 


U11748 


- 


+ 


+? 


irr, bright core/bar? 


U11819 





+ 


- 


fuzzy 


E0140040 


+ 


+ 


- 


spiral 


E0840411 


- 


- 


? 


edge-on 


E1200211 


+ 


- 


+ 


fuzzy mag. bar 


E1870510 


+ 


- 


+ 


irr. spiral, flocc. 
spiral 


E2060140 


+ 


+ 




E3020120 


+ 


+ 


+? 


spiral, hint of bar? 


E3050090 


+ 


+ 


+ 


barred spiral 


E4250180 


+ 


+ 


- 


spiral 


E488-049 


+ 




+ 


inclined mag. bar 



'^"-1-" indicates good fit p > 0.95; "0" indicates average 
fit 0.05 <p< 0.95; "-" indicates bad fit p < 0.05. 



^"+" indicates component clearly present; 
component not obviously present. 



indicates 



Table 11 
Morphology and probability NFW halos 



quality 



bar 



core 



both 



good {p > 0.95) 

bad {p < 0.05) 

unclear (0.05 <p < 0.95) 



4 


7 


3 


2 


4 





4 


7 


1 



Note. — Value indicates number of minimum 
disk fits of that quality in presence of component 
mentioned. 



